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ABSTRACT 

A sample of 427 gamma-ray bursts from a database (February 2002 - April 
2008) of the RHESSI satellite is analyzed statistically. The spectral lags and 
peak-count rates, which have been calculated for the first time in this paper, 



-2- 



are studied completing an earlier analysis of durations and hardness ratios. The 
analysis of the RHESSI database has already inferred the existence of a third 
group with intermediate duration, apart from the so-called short and long groups. 
First aim of this article is to discuss the properties of these intermediate-duration 
bursts in terms of peak-count rates and spectral lags. Second aim is to discuss 
the number of GRB groups using another statistical method and by employing 
the peak-count rates and spectral lags as well. The standard parametric (model- 
based clustering) and non-parametric (K-means clustering) statistical tests to- 
gether with the Kolmogorov-Smirnov and Anderson-Darling tests are used. Two 
new results are obtained: A. The intermediate-duration group has similar proper- 
ties to the group of short bursts. Intermediate and long groups appear to be dif- 
ferent. B. The intermediate-duration GRBs in the RHESSI and Swift databases 
seem to be represented by different phenomena. 

Subject headings: gamma-ray burst: general 



1. Introduction 



Mazets et al. I (Il98ll ) ; iNorris et al. I ( 1l984j ) ; iKouveliotou et al. I (Il993[ ) ; and lAptekar et al. 



( 1l998l ) have suggested the division of gamma-ray bursts (GRBs) into two categories, ei- 
ther short or long, according to their duration (at ~2s). Many observations demonstrate 
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D08[) . At present there is a predominant opinion that they are physically 
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There are also statistical indications of a third, "intermediate", group. The divi- 
sion of GRBs into thre e groups has been studied statistically oyer different data 



oases: 



BATSE fHorvath 



Horvath et al. 



2008 
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2001 



Horvath 



200d:IChattoDadhvav et al. IboOTh: B eppoS AX faorvath I l2009h: Swift JHorvath et al" 



2002 



Veres et al. 1120101 ) and RHE SSI flRipa et al. 



2009 ) . These three group s may also have different celestial distributions (IMeszaros et al. 



Vavrek et al. 2008[ ): at least for the BATSE database. No test has given a statistically 



significant support for the existence of four or more g roups. Only the BAT SE database gave 
a weak 6.2% significance level for such a possibility (iHorvath et al. 1120061 ). 
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One cannot e xclude an eventuality that the separation of this t hird group is simply 



a selection effect (IHakkila et al. I l2000l : iRajaniemi fc Mahonen II2002I ). In other words, a 



separation from the statistical point of view does not necessarily also indicate astrophysically 
different phenomena. In principle, it is still possible that the class of the intermediate 



GRBs constitutes a "tail" of either the short or the long group. The artic le by I Veres et al. 



( I2OIOI ) claims that - at least for the Swift database (ISakamoto et al. 1120081 ) - the third group 



is related to the so-called X-Ray Flashes (XRFs) , whi ch need not be physically distinct 



phenomena (IKippen et al. Il2003t ISoderberg et al. 



20061 ). Tw o models of XRFs are favor ed; 



either they are ordinary long GRBs viewed sli ghtly off-axis (Zharig fc Meszaros Il2002[ ) or 



they are intrinsically soft long-duration GRBs (iGendre et al 



20071 ). Hence, at least in the 



Swift database, the problem of the intermediate class seems to have been solved. 

However, for three reasons the situation has not yet been clarified. First, with regard 
to the Swi ft database, another study suggests that even the short group should be further 
separated (ISakamoto fc Gehrels II2009I ). Secondly, there is additional observational evidence 
against the simple scheme that maintains the existence of only two types of bursts (short/hard 
and long/soft) separated at duration of ~2s. The GRB060614 event, which is clearly long at 
duration (~ 100 s) but in any other properties resembles a short GRB, and subsequent short 



bursts with soft extended emission yet challenged this scheme (iGehrels et al. Il2006[ ). To 



avoid the limitations of a short -long separat i on terminology, the designations "Typ e I" and 



"Type H" have been proposed (jZhang 



20061 : IZhang et al. II2009I : iKann et al. Il2011f ) because 



duration alone is hardly sufficient for a correct division into categories. Thirdly, it remains 
possible that in other databases the discovered intermediate group is not represented by 
XRFs. Concerning this third reason, the mean duration of the intermediat e group appears 



to vary according to th e database in which it is found. For the Swift data (IHorvath et al. 
20081 : iHuja et al. 1 12009| ) the mean duration is ~ 12 s, which resembles the durations of the 



long GRBs, but for the RHESSI and B ATSE data flHorvath Ill998l : iMukheriee et al. Ill998 
Horvath et al. II2006I : iRipa et al. 1 120091 ) this mean is far below 10 s. 



It is clear that any new result in the classification scheme of GRB groups i s desir able. 
In this article we study the RHESSI database, where in addition to iRfpa et al. I ( 120091 ) . the 
spectral lags and peak-counts are also included. We have two concrete aims here: First, to 
provide further statistical tests concerning the GRB classes and, second, to provide additional 
inf ormation concernin g the physical significance of the RHESSI intermediate group found 
bv lRipa et al~l J2009h . 



The paper is organized as follows: In Sec. [2] the RHESSI satellite and its GRB data 
sample are described. In Sec. [3] distributions of spectral lags, normalized lags, and peak- 
count rates are studied using Kolmogorov-Smirnov and Anderson-Darling tests along with 
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Monte Carlo simulations. In Sec. H] we discuss results of these tests, compare the results 
with the BATSE and Swift data samples, and discuss the number of GRB groups using 
model-based and K-means clustering methods. Sec. E] summarizes results of this paper. 



2. The RHESSI data sample 



The Ramaty High Energy Solar Spectroscopic Imager[Jo (RHESSI ] is a satellite des igned 
for the observation of hard X-rays and gamma-rays from solar flares (ILin et al. II2002I ). but 
it is also able to detect GRBs. There is no automatic search routine for GRBs, and only if a 
message from any other instruments of the International Planetary Network (IPN) occurs, 
the RHESSI data are searched for a GRB signal. Therefore, our data set includes only events 
confirmed by other satellites. 



In this paper we study the same list of bursts which has been published in iRfpa et al. 



(]2009f ). We consider 427 GRBs from period between February 14, 2002 and April 25, 2008. 



Contrary to iRfpa et al. I (|2009[ ). the spectral lags and the peak-counts - calculated for the 
first time for RHESSI - next to the durations and hardnesses, are used. They are collected, 
together with their uncertainties, in Table |9l These new observational data allow further 
study of the questions concerning the GRB classifi cation. There are two arguments for 



the choice of the same list of bursts. First, both in iRfpa et al. I (|2009[ ) and in the present 
work, similar statistical studies are performed. Hence, for comparison, it is reasonable to 
study the structure of groups found in the RHESSI database over the same set. The second 
argument concerns an instrumental effect. The measurements of the hardness ratio of the 
events during the year 2008 and later has been systematic ally affected by an "annealing" 
pro cedur^ executed on the RHESSI detectors at late 2007 ( Bellm et al. 2008). The reason 
why the RHESSI team decided to anneal the detectors was to recover its deteriorating 
spectral sensitivity. However, the sensitivity at low energies had not been recovered as well 
as at high energies and hence the measured GRB hardness ratios from the post-annea ling 
period are systematically shifted to higher values (I Veres et al. II2009I : iRipa et al. Il2010[ ). In 
order to be eliminate this instrumental influence a more sophisticated modeling would be 
required. However, this is beyond the scope of this article. 



In order to compare the spectral lags and the peak-counts of bursts belonging to the 
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different groups, one must provide a rule by which the particular GRBs are sorted into the 
concrete groups. We proceeded in the following manner: The probabi li ty de nsity function 



employed in the fitting of the duration-hardness plane in iRfpa et al. I (|2009[ ) is composed 
through the summation of three bivariate log-normal functions f{x, y) = fi{x, y) + f2{x, y) + 
f3{x,y), where x is the base 10 logarithm of the duration and y is the base 10 logarithm of 
the hardness ratio, and /i, /2, and /s are components corresponding to the particular groups. 
A burst at the point [xq; yo] is considered short, intermediate or long depending on whether 
the f^(xo,yo), f2(x(),y()) or .f-j(xn,j7n) is m aximal. In essence, we follow a procedure identical 



(120061 ) and lHorvath et al. I (|2010|) utilized on the BATSE and 



to the that of iHorvath et al. 
Swift datasets. 

In order to sort the given GRBs into the grou ps we employ the m easurements of the 
durations and hardness ratios as given in Table 7 of iRipa et al. I ( 120091 ) with the exception 
of six events. We found that for these six events the mentioned values in iRipa et al. I (120091 ) 
were not corrected for a so-called decimation, which is an instrumental mode used to conserve 
the onboard memory. Table [1] presents these six events, now corrected for this decimation. 

The group r nembers used in th is study were determined from the best Maximum Like- 
lihood (ML) fit ( IRfpa et al. 1 12009[ ) in the duration-hardness plane of 427 GRBs. In this 
sample the six eve nts with corr e cted d ecimation were included along with the remaining 421 
events taken from iRipa et al. I ( 120091 ). The best ML fit with two bivariate lognormal com- 
ponents gives logarithmic likelihood lnL2 = —313.4. Best fit with three components gives 
logarithmic likelihood InL^ = —303.4. The ML ratio test tells that the twice of difference in 
the logarithms of t he likelihoods 2('lnL3 — lnL2) = 20.0 should follow distrib ution with 6 
degree s of freedom (IHorvath et al. 1120061 ). Therefore ML ratio test, employed in lRipa et al. 
( 120091 ) and now applied on the duration-hardness plane with these 427 GRBs including the 
six events corrected for decimation, gives again a statistically significant intermediate group 
at the significance level of 0.3 %. The new (former) best-fit model parameters of the inter- 
mediate group are: 0.12 (0.11) for the mean logarithmic duration, 0.25 (0.27) for the mean 
logarithmic hardness, 4.1 (5.3 %) for the weight, and 0.0 (0.59) for the correlation coefficient. 
The group members are shown in Fig. [1] and listed in Table |9l 

The spectral lags L of the RHESSI data were calculated by ourselves by fitting the peak 
of the cross-correlation function (CCF) of the background-subtracted count light-curves at 
two channels, 400 — 1500 keV and 25 — 120 keV, by a third order polynomial. The position 
of the maximum of the polynomial fit measures the spectral lag. An example of such a 
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Fig. 1. — The hardness ratio H plotted against the duration Tgo for the RHESSI database 
with the best ML fit of three bivariate log-normal functions. The different GRB group mem- 
bers are denoted with different symbols: the crosses, full circles, and triangles correspond, 
respectively, to the short, intermediate, and long bursts. CL means "confidence level". 

calculated for the RHESSI GRBs. 

To obtain statistical errors, a Monte Carlo (MC) method was utilized. The following 
procedure was employed to prepare 1001 synthetic count profiles for each GRB: The measured 
count profiles were randomly influenced by Poisson noise, after which the background was 
subtracted. The RHESSI count rates are sometimes "decimated", which means that, as the 
rate becomes too high or the onboard sohd-state recorder becomes too full, a part of the 
recorded counts is removed. If decimation occurs, the fraction (/d — 1)//^ of the counts below 
a decimation energy Eq is removed, /^is the decimation factor (weight), usually equal to 4 or 
6. All events above Eq are downlinkeco. To prepare the synthetic count profiles, the number 
of counts in each bin was changed accordingly to the Poisson distribution. The 1-sigma 
errors for non-decimated, fully decimated, and partially decimated data are y/C, -y/Zd-Cdc, 
and a/Ci + /d-C*2,dc) respectively. C is the measured count number in a bin for non-decimated 
data. Cdc is the count number in a bin of fully decimated data and consequently corrected 
for this decimation. Ci is the count number in the non-decimated portion and C2,dc is the 
corrected count number in the decimated portion of the measured rate in the case of partially 
decimated data. A detailed explanation is provided in Appendix |X1 The CCF was fitted for 
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Fig. 2. — Left: An example of the cross-correlation function of two background-subtracted 
count light-curves of the very bright GRB 060306 derived at two energy bands 400 — 1500 keV 
and 25 — 120keV. Right: A detail of the same curve with the third order polynomial fit (thick 
solid curve). The position of the maximum of the fit measures the spectral lag (dotted line). 
The boundaries of the polynomial fit are marked with dashed lines. 

each of the 1001 synthetic profiles and for each burst in our sample. The median of such a 
distribution of 1001 maxima of polynomial fits was taken as the true lag L for each burst. 
Theses median lags L are used in the following statistical tests and listed in Table [H The 
2.5% and 97.5% quantiles of such a distribution of 1001 maxima of polynomial fits for each 
GRB delimit the 95 % CL statistical errors. These errors are also listed in Table M 

We decided to calculate the spectral lags only for bursts with a signal-to-noise ratio 
higher than 3.5 in both channels. This signal-to-noise ratio is defined as ^Tgo/ V'S'tqc + 2-Bt90) 
where 5x90 is a GRB signal over the background level -BtqO) and both S and B are counts 
in a Tgo time interval over the range 25 keV— 1.5 MeV. The choice of this limit was found 
to ensure that the CCF was sufficiently smooth with a clear peak allowing determination of 
a reliable lag. Therefore, excluding the noisiest data, the number of GRBs with calculated 
lags is 142. Their distribution is presented in Fig. |3l 

The GRB peak-count number 5* was derived from the light-curve with the maximal 
count number C at the range 25keV— 1.5MeV after subtracting the background B. The 
peak-count rate F is given as the peak-count number S divided by the width of the time bin 
Stj-es- This width was different for different GRBs, and covered a range between 2 ms and 
3 s. The dimension of the peak-count rate is count / s. 

The one sigma error a-p of the peak-count rate F was calculated as a-p = cs/Stj-es, where 
the error of the GRB peak-count number is as = a/ (o"c)^ + (c"b)^- We assume that errors 
of the maximal count numbers cxc and of the background Ub = \/B are Poissonian and 
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independent. The error ac is: crc = y/C in case of non-decimated data; given by expression 
flA2p in case of fully decimated data; and given by expression flA4p in case of partially 
decimated data (see Appendix The peak-counts with errors were calculated for all 427 
objects. 



3. Properties of the GRB groups 



3.1. Distribution of spectral lags 



In this section we use Anderson-Darling (A-D) test (lAnderson fc Darling Ill952l : iDarling 



19571 ) to compare distributions of spectral lags of different GRB groups (see Fig. [3]) found 
by the ML method applied on durations and hardness ratios (see Sec. [2]). The short (inter- 
mediate, long) group contains 26 (11, 105) objects. The mean values of the spectral lags of 
these groups are similar, hence we use the A-D test becau se it is particularly sensitive to the 
tails of the tested distributi ons (IScholz fc Stephens 1ll987l). For its calculation we employ adk 



201ll ). The results are summarized 



package of the R softwar^ ( R Development Core Team 
in Tab. [3 

The A-D test gives a significance of 16.8% (the probability that the two samples are 
drawn from the same distribution) for the short-intermediate pair, and it yields a significance 
of 4.2 % for the long-intermediate pair. Therefore, in case of short and intermediate groups, 
we cannot reject the null hypothesis that the two samples are drawn from the same distribu- 
tion on a sufficiently low level (5 %). On the other hand, this null hypothesis can be rejected 
in the case of long and intermediate groups, but the significancy is not far below 5 % level. 
The same test applied on the lags of the short-long pair yields a significance of < 10"'^ %. 
Therefore, in this case, the null hypothesis can be rejected with a high significance. This 
strongly supports the well-known claim that the shor t and long GRBs are really different 
phenomena and confirms the results of iNorris et al. I ( 120011 ) (obtained with BATSE), but 
now by using the RHESSI instrument. 



3.2. Distribution of normalized lags 



In this section we compare the distributions of normalized lags (Fig. H]), i.e. L/Tqq, 
next to the absolute values of the lags. Again we use the A-D test between the different 
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GRB groups mentioned in the previous section. The number of events within the groups is 
therefore the same. The results are summarized in Tab. [31 

The A-D test gives the significance level of 54.2 % for the short-intermediate pair and 
it gives the significance level of 45.0 % for the long-intermediate pair. The significances are 
considerably above 5 % level, therefore the null hypothesis that the samples are drawn from 
the same distribution cannot be rejected. For the short-long, pair the A-D test gives the 
significance level of 6.0 %. If the normalized lags are concerned, the difference between the 
short and long bursts is not definite. 



3.3. Distribution of peak-counts 



Here we used Kolmogorov-Smirnov (K-S) test ( Kolmogorov 1933 : Smirnovlll948l ) to 
compare the cumulative distributions of the peak-counts among the different GRB groups. 
The short (intermediate, long) group contains 42 (18, 367) objects. The results are presented 
in Tab. [Hand shown in Fig. [51 

The results of the K-S tests imply that the distributions of the peak-count rates are 
different over all three groups. Particularly, the K-S significance level for the intermediate 
vs. short busts is 0.9%, for intermediate vs. long bursts it is 3xl0~^%, and for short vs. 
long bursts it is < 10~^ %. 



3.4. Monte Carlo simulations 

In order to test the robustness of the results obtained by the A-D tests applied on lags 
L, normalized lags L/Tgo and K-S tests applied on peak-count rates F, one can use Monte 
Carlo method. 

In case of spectral lags we proceeded in the following way: The procedure described in 
Sec. [21 - calculation of statistical errors of the lags by applying of Poisson noise - provided 
distribution of 1001 lags for each GRB. Thus for each GRB we randomly selected one lag 
from its distribution and made 10 000 data samples. Then the A-D tests for these 10 000 
samples were calculated. 

In case of peak rates, we proceeded as below: We applied the Poisson noise to the 
measured light curves and subtracted the background in order to obtain the simulated data. 
Then we derived the peak count rate for the same peak time when the peak was found in 
the measured light curves. We proceed in this way for each GRB. Afterwards we calculated 
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K-S tests and repeated this sequence 10 000 times. 

The number of cases when the A-D and K-S probabihty reached higher values than 
5 % for tests done on different pairs of GRB groups is noted in Tab. O The results of MC 
simulations comparing spectral lags and normalized lags are shown in Fig. [61 

The MC method confirms that the distributions of spectral lags between short and 
long GRB groups are different. Let's compare results from the MC simulations of lags and 
normalized lags between the intermediate-short and intermediate-long pairs with the results 
of the tests applied directly on median lags (Tab. [2]) and median normalized lags (Tab. |3]). 
Then one can see that MC simulations gives A-D prob. > 5 % more often then expected. 
It can be caused by the fact that for some GRBs, weak and noisy ones, the distribution of 
lags found by MC method might not follow the real distribution because after applying the 
Poisson noise the polynomial fit of CCF may not well describe the CCF peak. The reason for 
this conjecture is that the fitting range remained fixed and same for the simulated data as 
for the measured data. In other words the suitable fitting range for the measured data need 
not be suitable for the simulated data. In this case we think that the A-D tests applied on 
the median lags give more reliable results than the MC simulations do. However, one mutual 
behaviour is seen here, the intermediate-short pair has distributions of lags and normalized 
lags more similar than the pair intermediate-long does. This feature is seen both in the 
A-D tests applied on the median lags/normalized lags and in the A-D tests of the MC data 
samples. 

MC simulations also confirm results of K-S tests applied directly on the measured peak 
rates. We can conclude that the short-, intermediate- and long-duration bursts have different 
distributions of peak count rates. The results of MC simulations comparing spectral lags 
and normalized lags are shown in Fig. [71 



4. Discussion 



4.1. Comparison with the BATSE database 



The lags of GR Bs from the BATSE dataset are different for the short and long groups 
(iNorris et al. 1120011 ): for the short bursts the lags on average are close to zero, but for the 
long bursts they are positive. Norris and his collaborators did not study the lags of th e 



intermediate bursts separately ( iNorris et al. I l200ll : iNorris I l2002l : [Norris fc Bonnell II2006I ). 



For the sake of compl eteness we have attempted to do this for the publicly available data. 



Horvath et al. I ( |2006[ ) defines membership within the groups for all BATSE GRBs. Addi- 
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tionally, iNorris I ( 20021 ) defines for any GRB witli T90 > 2 s its lagj. Compilation of tliese 
two lists and tlie application of the A-D test on the lags of the three BATSE groups (here 
for the first time) produced the results collected in Tab. M and shown in Fig. |H1 The short 
(intermediate, long) group contains 33 (119, 1179) objects here. Of course, one must keep 
in mind that this sample is drastically truncated for the short bursts. Hence, the short- 
intermediate and the short-long comparisons, can serve only as a qualitative indicators. 
Even the intermediate-long pair cannot be taken as representative because the truncation 
T90 > 2 s can also omit several intermediate GRBs. 

Keeping all this in mind, if the lags are taken into consideration, we can say that there 
is some similarity between the BATSE and the RHESSI databases. First of all, there is 
a similarity with regard to the intermediate-long pair: the difference is confirmed, though 
not on an high significance level, but remarkably the significances from the A-D test are 
comparable (3.8% and 4.2%). Second, there is a similarity with regard to the intermediate- 
short pair: in both databases the A-D test reveals that for these two groups the distributions 
of GRB lags are similar; the significances are 51.3% for BATSE, and 16.8% for RHESSI. 
However, one must again keep in mind that our BATSE sample of short bursts is truncated 
and the sample of intermediate bursts as well. Third, both databases show a difference 
between the average lags for the short-long pairs (for the BATSE databases the difference 
between the distributions is shown not to be signi ficant; the A-D P-y alue is only 9.7%, 
probably as a result of the sample truncation, but iNorris et al. I (120011 ) makes this claim 
unambiguously) 

The results of the K-S tests applied on the peak-fiuxes of the 64 ms resolution light- 
curves for the BATSE data imply that the distributions are different over all three groups. 
The short (intermediate, long) group contains 502 (169, 1282) objects here, and the K-S 
tests are sum med in Tab. [7| and shown in Fig. [HI These results are to be expected because, 
for example, iNakar I (120071 ) claims that the peak-fiuxes of short GRBs are roughly 20 x 
smaller than those of the long ones. It is also known that the i ntermediate BATSE group is 
"intermediate" concerning the fiuence (IMukherjee et al. Ill998l ). 



Therefore, our comparison of these RHESSI and BATSE groups finds similarities. In the 
case of BATSE database, all three groups are different in respect to two quantities (duration 
and peak- flux). It is remarkable that fo r BATSE the hardne ss of intermediate group is 
strongly anticorrelated with the duration (IHorvath et al. II2006I ). Since the hardness of the 
intermediate group differs from the hardnesses of the short and long ones, these studies 
support the opinion that all three BATSE groups represent different phenomena. 



^l\protec-t||ittp://heas arc. gsfc.nasa.gov/docs/cgro/analysis/lags/web Jags. html 
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4.2. Comparison with the Swift database 



The lags of the GRBs from the Swift dataset are al so different for the short and long 



groups (Ide Ugarte et al. Il201ll ). Ide Ugarte et al. I (120 111 ) also discuss the lags of the inter- 



mediate bursts, and they find a behavior which does not resemble the cases found in the 
RHESSl and BATSE datasets. The Swift's intermediate-long pair has on average similar 
lags, but there is a statistically significant difference in the short-intermediate pair. Thus, if 
the lags are considered , the Swift's intermediate group is similar to the Swift's long group 
(Ide Ugarte et al. 11201 ll ). On the other hand, the peak- fluxes differ signiflcantly in the short- 
intermediate and intermediate-long pairs, respectively. The peak-fluxes o f the short-long 
pair are not different f rom the statistical point of view ( IVeres et al. Il2010l ). Nevertheless, 
de Ugarte et al. I (120 111 ) concludes that "Swift's intermediate bursts differ from short bursts, 
but exhibit no signiflcant differences from long bursts apart from their lower brightness" . In 
other words, in the Swift database there is a clear similarity between the intermediate group 
and the long o ne. The physical difference of the short an d long bursts in the Swift database 



further holds (IVeres et al. II2OIOI : Ide Ugarte et al. 1120111 ) 



Comparison with the Swift's groups leads to the conclusion th at the third g r oup i n the 



Swift database is stron gly related to the long group, as stated by I Veres et al. I (I2OIOI ) and 



de Ugarte et al. I ( 1201 it ), and only the short group should represent another phenomenon. 



There is a difference in the hardness, peak- flux, and duration for the intermediate-long 



pair (iHorvath et al. 1 120081 : IVeres et al. Il2010[ ). but no clear separation occurs for the lags 



(Ide Ugarte et al. Il201l[ ). We have no reason to query the conclusions of I Veres et al. I (I2OIOI ) 



that the intermediate group is related to XRFs which in turn can be related to standard 
long GRBs. We add th at the separation within the long grou p itself into harder and softer 



parts is not fully new (jPendleton et al. I Il997l : iTavani Ill998l ). We allow to claim that the 



intermediate-duration bursts in the RHESSl and Swift databases are different phenomena. 
These results followed exclusively from the statistical analyses. 



4.3. Discussion of the number of groups 



In order to provide an extended discussion on the number of GRB groups we apply 
clustering met hods to our data sample. This also serves to extend the statistical analysis 
performed by iRipa et al. I (120091 ). In general, the clustering methods can be divided into 
parametric and non-parametric types. Parametric methods assume that the data follow a 
pre-deflned model (in our case a sum of multivariate Gaussian functions). These methods 
assign for each GRB a probability of membership in a certain group. The non-parametric 
methods, e.g. K-means clustering, provide deflnite assignments of each burst to a given 
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group. More details about these methods c an be found in the book by lEveritt et al. 
Model-based clustering is also described in iMcLachlan &: Peel I ( 120001 ) . 



(1201 ih . 



We apply model-based clustering and K-means clustering methods on our RHESSI data 
sample by using the algorithms implemented in the R software. 



4-3.1. Model-based clustering method 

In this method we assume that the distribution of the tested parameters (logarithms of 
durations, hardness ratios, peak-count rates, and normalized lags) follow a superposition of 



(2006 


) and 


Veres et al. 


(2010) 



The Maximum Likelihood method is used to find the best-fitted model parameters. 
Adding more free parameters to a fitted model can increase the likelihood, but also may 
result in overfitting. It is possible to penalize a model for more free parameters. T his can be 



done by a method called the Bayesian Information Criterion (BIG) presented by ISchwarz 
( 119781 ). The function which must be maximized to get the best-fitted model parameters 
is: BIG = 21nlmax - mlnN, where Imax is the maximum likelihood of the model, m is the 
number of free parameters, and N is the size of the sample. In our work we use the BIG to 
determine the most probable model, its parameters and the number of its components. 

For model-based clustering, we use Mclust packag (0 JPralev fc Rafterv I boooh of R. 
For the explanation of different models, see the Mclust manuaj^. The nomenclature of the 
different models in Mclust involves the following designations: the volumes, the shapes and 
the orientation of the axes of all clusters may be equivalent (E) or may vary (V) and the 
axes of all clusters may be restricted to parallel orientations with the coordinate axes (I). 



4.3.2. Model-based clustering - 2 variables 
First, we start with a two-dimensional case and fit Tgn durations and hardnesses H. The 



data sample consists of 427 bursts (Table 7 of iRipa et al. I (120091 ) and Tabled]). 

In this case the number of free parameters of the model with k bivariate Gaussian 



^http://cran.r-project.org/web/packages/niclust/index.html 
^http://www. stat.washington.edu/research/reports/2006/tr504.pdf 
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components is 6A; — 1 {2k means, 2k standard deviations, k correlation coefficients, and k — 1 
weights, because the sum of the weights is 1). For the most general model, all parameters 
are free. However, sometimes we want to test models where some of the parameters between 
different components are in a relation with other parameters, e.g. all components have the 
same weight or shape, etc. In this case the number of degrees of freedom is reduced. 

As seen in Fig. |9] the best fitted model has k = 2 components with equal volumes, 
variable shapes, and with the axes of all clusters parallel with the coordinate axes (EVI 
model). This best fitted model has a value of BIG = —681.5. The EVI model with one 
component gives BIG = —899.1 and with three components BIG = —701.8. For all other 
tested models with k = 1 component the highest BIG is -820.3 and with k = 3 components 
-694.3, which are clearly below the maximum. 

The diff erence between the BIG of t wo models gives us inform ation about the goodness. 



According to lKass fc Raftery I (119951 ) and lMukherjee et al. I ( 119981 ) . a difference in BIG of less 



than 2 represents weak evidence, difference between 2 and 6 represents positive evidence, 
between 6 and 10 strong evidence, and difference greater than 10 represents very strong 
evidence in favor of the model with the higher BIG. 

In our case the difference between the best fitted model (EVI) with two components 
and the EVI models with one or three components is always higher than 10. This gives a 
strong support for the EVI model with k = 2 components. 

The two components are short/hard and long/soft groups. The intermediate-duration 
bursts showed in Fig. [T]are assigned to the short /hard group by this test. 



4-3.3. Model-based clustering - 3 variables 

Next we perform model-based clustering on three variables: Tgg durations, hardnesses 
H, and peak-count rates F. Since the peak rates were measured for all events, the sample 
here also consists of all 427 bursts (Table |9]). 

The best fitted model has k = 3 components (see Fig. [9]) with equal volumes, equal 
shapes and equal correlation coefficients between all clusters (EEE model). This best model 
has a value of BIG = —1156.6. The EEE model with two components gives BIG = —1168.7 
and for four components BIG = —1174.6. Markedly high values of BIG are also obtained for 
the EEI, VEI, and VVI models with k = 3 components, BIG = -1166.2, BIG = -1162.5, 
and BIG = —1166.1, respectively. 

The difference in BIG between the EEE model with three components and the EEE 
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models with two or four components is > 10. The other models with other number of 
components (except above-mentioned EEl, VEl, and VVI models with three components) 
gives BIG value lower by at least 10. This provides a strong evidence in the favor of EEE 
model with k = 3 components. The group structure of this model with three components is 
shown in Fig. [TOl 

The intermediate-duration bursts showed in Fig. [1] are assigned to the short/hard group 
by this test. A new result here is that the group of long bursts is separated into high- and 
low-peak flux clusters. 



4.3.4- Model-based clustering - 4 variables 



Spectral lags of BATSE GRBs, i.e. the time delay between low and high energy photons 
from short and long groups have been found to differ. For short bursts, an average lag is 
~ 20 — 40 times shorter than for long bursts, and the lag distribution is close to syrnmetri c 
about zero - unlike long bursts ( iNorris et al. 1 1200 ll : iNorris 1 120021 : iNorris fc Bonnell 1 12006| ) . 
This result gave us an idea to incorporate the spectral lags as well. 

In this part we apply the model-based clustering on GRB peak-count rates F, Tgo 
durations, hardness ratios H, and as a new addition to the variables, on normalized lags 
L/Tgo. Since the RHESSI spectral lags were calculated only for 142 bursts (Table [9]), our 
sample is truncated. 

The best fitted model has k = 2 components and is unconstrained, i.e. it has variable 
volumes, variable shapes and variable correlation coefficients (VVV). The best BIG value 
for this model is BIG = —1768.4. However, the VVV model with k = 3 components gives 
similar value of BIG = —1768.5. The other models gives BIG value lower by at least 10. This 
strongly supports the VVV model with k = 2 components. There is no need to introduce 
VVV model with three components that has more free parameters. The two components are 
separated accordingly to the values of normalized lags into zero- and non-zero-lag events. 



4.3.5. Summary of model-based clustering 

The model-based clustering of two-parameter data (Tgo and H) gives a strong evidence 
in favor of EVI model with two components. The analysis of three-parameter data (Tgo, 
H and F) shows that the best fitted model is EEE with three components. Surprisingly 
a new result is obtained here; the group of long bursts is separated into high- and low- 
peak fiux clusters. The analysis of four-parameter data (Tgo, H, F and L/Tgo) supports the 
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VVV model with two components only. The separation into the two components here is 
accordingly to the values of normalized lags into zero- and non-zero-lag events. 



4-3.6. K-means clustering 



One of the non-parametric clustering methods is K-means (IMacQueen 1119671 ). Before we 
use our data for this method we scale them, i.e. we subtract the mean value and then divide 
them by the standard deviation. The reason for this procedure is that the clustering algo- 
rithm is sensitive to the distance scale of the variables. For more detail s about the application 



of the K- means metho d in a similar analysis of GRB data, see, e.g. IChattopadhyay et al. 



( 120071 ) or IVeres et al. I (120101 ). For this clustering method we use kmeans package imple- 
mented in the R software. 

To use the K-means method, one must set the number of clusters beforehand. Then 
the corresponding number of centers is found by minimizing the sum of squared distances 
from each burst to the center of the group to which they belong. There is no precise way 
to determine the best number of clusters with this method. However, it has been suggested 
that if one plots the within-group sum of squares (WS S) as a functiori of the number of 
clusters, then an "elbow" will indicate the best number ( IHartigan Ill975l ). This method do 
not provide any probability indicating the significance or insignificance for the given best 
number of clusters. 

The calculated WSS as a function of the number of groups for our data samples using 
2 (Tgo, H), 3 (Tgo, H, F), and 4 (Tgo, H, F, L/Tqq) variables are rather smooth and do not 
demonstrate any remarkable and sharp "elbows" and thus do not bring useful information 
on the GRB classification. 



4.4. Discussion of the results 

The K-S tests applied on peak-count rates show that the distributions are different over 
all three groups. The K-S significance level for the short-long pair is < 10^^ %, for the 
intermediate-long pair it is 3x10^^%, and for the intermediate-short one it is 0.9%. The 
short and long GRBs have clearly different distributions of peak rates. Also the intermediate 
and long GRBs have clearly different distributions of peak rates. The intermediate-short pair 
also exhibits different distributions (K-S probability < 5 %), however less markedly than the 
other pairs of groups do. These results are confirmed by MC simulations. 

The A-D tests applied on distributions of spectral lags unveil that the A-D probability 
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for the short-long pair is < 10~^ %, for the intermediate-long pair the A-D probability is 
4.2 %, and for the intermediate-short one it is 16.8 %. The short and long GRBs have 
clearly different distributions of spectral lags. The intermediate and long GRBs have A-D 
probability < 5%, however in this case the difference is not strong. The intermediate-short 
pair does not exhibit different distributions. The difference in the spectral lag distributions of 
the short-long pair of ORB groups is confirmed by MC simulations. In case of intermediate- 
short and intermediate-long pairs the MC simulations reveal the same tendency as the A-D 
tests applied directly on the measured values, i.e. the intermediate-short pair has more 
similar distributions of spectral lags than the intermediate-long pair has. However, MC 
simulations gives A-D probability higher than 5 % more often then expected. A possible 
reason is commented in Sec. 13.41 

The A-D tests applied on distributions of normalized lags show that these distributions 
cannot be claimed as different. The A-D probability for the short-long pair is 6.0%, for the 
intermediate-long pair it is 45.0%, and for the intermediate-short one it is 54.2%. Here one 
can see the same tendency as in the case of A-D tests applied on spectral lags, i.e. the short- 
long pair has the least similar distributions of lags, the intermediate-long couple stays in the 
middle, and the intermediate-short pair has the most similar distributions. This tendency 
appears also in the MC simulations, however the absolute frequency of the cases when A-D 
probability exceeds 5 % level happens more often then expected. The reason could be the 
same as in the case of MC simulations applied on absolute values of spectral lags. 

The model-based clustering of two-parameter data (Tgo and H) gives a strong evidence 
in favor of a model with two components only. The two components are short/hard and 
long/soft groups. The intermediate-duration bursts showed in Fig. [1] are assigned to the 
short/hard group by this test. The analysis of three-parameter data (Tgo, H and F) shows 
that the best fitted model has three components. Surprising point here is that this method 
separates the group of long bursts into high- and low-peak flux clusters. The analysis of 
four-parameter data (Tgo, H, F and L/Tgo) supports a model with two components which 
are separated accordingly to the values of normalized lags into zero- and non- zero-lag events. 

Surveying the A-D and K-S tests of the RHESSI data, it should be raised that the 
difference between the short and long bursts was again strongly confirmed. This follows 
from the different distributions of the spectral lags and by the different distributions of 
the peak-count rates, both results were confirmed by the MC method. This is already an 
expectable result, but - usefully - this result came from a new observational database. 

According to Fig. [T] the intermediate-short pair of groups have similar hardness ratios. 
Also according to the results of the A-D test of the spectral lags, the distributions of lags are 
not different for the intermediate-short pair. However, the intermediate-duration and short- 
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duration bursts are not completely same because their peak-count rate distributions differ. 
On the other hand, the intermediate-long pair of groups differs in hardness ratios, spec, lags, 
and peak-count rates. Therefore, in ou r opinion, it is po ssible that the intermediate group 
detected by RHESSI in Sec. [2] and by (IRipa et al. 1 120091 ) may be a longer tail of standard 
short/hard bursts. This can be supported also by the fact that model-based clustering 
method applied on hardness ratios and durations unveil only two clusters as the best solution; 
classical short/hard and long/soft groups, the intermediate-duration bursts are assigned to 
the short group. 

The RHESSI intermediate and long groups seem to be different phenomena. This differ- 
ence is supported by the distribution of the peak-count rates and spectral lags. The results 
show that the intermediate group is also "intermediate" with regar d to its lags. The inter- 
mediate group detected by Swift was found to be relat ed to XRFs (IVeres et al. Il2010l ). and 
those may in turn belong to the standard long GRBs ( Kippen et al. 20031 ). In the case of 
RHESSI, the longer and softer GRBs are more difficult to be detected, because RHESSI's 
sensitivit y declines rapidly below ^ 50 keV and the weak and soft GRBs are not easily ob- 
servable (IRfpa et al. 1120091 ). On the other hand. Swift is less sensitive in the photon-energy 
range > 150 keV. But softer GRBs are readily detectable with this instrument. Hence, in 
our opinion, an instrumental effect may be responsible for that the two satellites (Swift and 
RHESSI) detected different intermediate groups. This means that - from the statistical point 
of view - different groups can be found if one looks at different databases. 

There are bursts observed with properties simila r to the short bursts ( hard ness, lag) ex- 



cept t heir durations exceed 2 s. For exam ple, bevond iGehrels et al. I (120061 ) and lKann et al. 



( I2OIII ) mentioned already in Introduction. iNorris fc Bonnell I ( l2006l ) claims that "s hort bursts 
with extended emission (SGRBEE) can have Tqq > 2 s". Furthermore, also others (Ide Barros et al. 
20 111 ) propose the astrophysical fragmentation of the short GRB group. 



Concerning SGRBEE, we inspected the light-curves of all 18 RHESSI intermediate 
bursts, but we found no softer extended e mission comi n g afte r the main hard spike as is 
typical for this kind of bursts. Figure 3 of iPerley et al.l ( 120091 ) shows that the average Tqq 
duration of the initial spike of an SGRBEE lies between the average durations of short and 
long bursts. If RHESSI detects only the hard initial spike, and the softer extended emission 
is lost in the noise, then the detected intermediate group might be polluted by these objects. 
Therefore, we als o checked the light-c urves of the RHESSI intermediate bursts as observed 
by Konus-Wind (lAptekar et al.lll995l ). because it also has a good sensitivity below 50keV 
(its range is 10 — lOOOOkeV). It has an overlap with the following RHESSI intermediate 
bursts: GRB 020819A, GRB 030410, GRB 040329, GRB 050530, GRB 070802, GRB 070824, 
GRB 080408. However, no extended emission was observed for these seven bursts by Konus- 
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Wind. This observation indicates the RHESSI intermediate GRBs should not be dominantly 
polluted by SGRBEEs. 

Furthermore, there are also additional indications that GRBs which do not belong to the 
long+XR F pair category, rn a y orig inate from a broad range of astrophysical phenomena. For 
example, iMukherjee et al. I (119981 ) found four subclasses in the BATSE database from the 
year 1998, but the fourth group was populated by a single GRB. From the statistical point of 
view, such an object is an outlie r of uncertain origin. Likewise, similar sit uations exist con 



cerning the objects GRB 060614 flGehrels et al. 1120061 ) and GRB 110328A (ICummings et al. 



201ll ). Any study of such a single unusual object is beyond the scope of this article providing 
only statistical analyses. 



5. Conclusions 

The main results of this study can be summarized as follows: 

» Maxi mum Likelih o od te st in the duration-hardness plane of 427 RHESSI GRBs, taken 



from iRipa et al. I (|2009[ ) but now with six events corrected for decimation, again ex- 



hibits st atistically significan t third, in duration intermediate, group. This completes 



the work iRipa et al. I (|2009[ ) using the durations and hardnesses only. 



The spectral lags and peak-count rates have been calculated for GRBs observed by the 
RHESSI satellite for the first time. The spectral lags were obtained for 142 objects, 
and the peak-counts for all 427 GRBs. Hence, we constructed a new observational 
database for this satellite. Then the three GRB subgroups were analyzed statistically 
with respect to these new spectral lags and peak-count rates. 

The difference between short and long groups has been confirmed. Usefully this result 
came from a new observational database. 

Kolmogorov-Smirnov and Anderson-Darling tests applied on spectral lags and peak- 
count rates indicate that the intermediate group in the RHESSI database might be a 
longer tail of the short group or at least has some common properties with this short 
group. Contrary to this, the intermediate and the long groups are different. 

The group of RHESSI intermediate-duration GRBs is not dominantly populated by 
SGRBEEs. 

The intermediate-duration bursts found in the RHESSI and Swift databases seem to 
be represented by different phenomena. 
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Table 1. Six RHESSI GRBs with corrected T90 durations and hardness ratios. 



GRB'' 


Peak time*^ 


T90 (s)'= 


Hardness ratio log H'^ 


030518B 


03:12:23.050 


{1.86±0.07)E+1 


(2.90±0.27)E-1 


030519A 


09:32:22.500 


{3.20±0.27)E+0 


{5.31±0.61)E-1 


031024 


09:24:14.350 


(4.30±0.17)E+0 


-(2.06±0.31)E-1 


040220 


00:55:15.800 


{1.80±0.07)E+1 


{9.39±2.72)E-2 


050216 


07:26:34.275 


(4.50±0.56)E-1 


{2.33±0.48)E-1 


050530 


04:44:44.900 


(2.40±0.26)E+0 


{2.41±0.63)E-1 



''RHESSI GRB number. 

''Peak time of the count light-curve in UTC. 

"^The uncertaint ies were calculated though the same procedure used in 
iRipa et al. 

'^The hardness ratio was defined as the ratio of GRB counts at two 
different bands, H = S{i20-i500)kev/'5{25-i20)keV- 



Groups 


A-D P 


Group 


Mean L 


Median L 


cr 




(%) 




(ms) 


(ms) 


(ms) 


Inter.-Short 


16.8 


Short 


4.9 


1.9 


16.7 


Inter.-Long 


4.2 


Inter. 


28.7 


5.9 


78.4 


Short-Long 


< 10-3 


Long 


178.0 


50.8 


874.9 



Table 2: Left part: Results from the A-D tests of the spectral lag distributions for the 
RHESSI database, are presented. The null hypothesis is that the two samples are drawn 
from the same distribution. P denotes the P- value of the test. Right part: The means, 
medians and standard deviations a of the lags are listed. 



Groups 


A-D P 


Group 


Mean 


Median 


a 




(%) 




L(ms)/T9o(s) 


L(ms)/r9o(s) 




Inter.-Short 


54.2 


Short 


21.3 


15.8 


63.3 


Inter.-Long 


45.0 


Inter. 


17.6 


5.5 


63.0 


Short-Long 


6.0 


Long 


10.2 


3.4 


32.0 



Table 3: Results of the A-D tests of the equality of the normalized lag distributions between 
different RHESSI GRB groups are listed. P denotes the P-value of the test. Right part: The 
means, medians and standard deviations a of the normalized lags are also mentioned. 
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Fig. 3. — Left panels: The spectral lags of RHESSI GRBs sorted along the y-axis with 

respect to the value of (+error + |-crror|) for short-, intermediate- and long-duration bursts. 
The median lags, for each GRB, were taken from the lags of 1001 synthetic background- 
subtracted count time profiles obtained by Monte Carlo simulations of the measured profiles 
that were randomly infiuenccd by the Poissonian noise. The error bars are composed of 
the 95 % CL statistical error and the profile time resolution. A positive lag means that the 
low-energy counts are delayed. Right panels: The cumulative distributions of the obtained 
median lags for the three groups of bursts are shown. 
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Fig. 4. — The cumulative distributions of the normahzed lags for the three RHESSI GRB 
groups. 
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Fig. 5. — Left panel: Peak-count rates F of RHESSI GRBs as a function of Tgo durations 
for the three GRB groups, identified by the analysis of the hardnesses and durations, are 
displayed. Right panel: Cumulative distributions of these peak-count rates F for the short-, 
intermediate-, and long-duration bursts are shown. 



Groups 


D 


K-S P 


Group 


Mean 


Median 


a 






(%) 




F(s-i) 


F(s-^) 




Inter. -Short 


0.44 


0.9 


Short 


9 485 


5163 


20418 


Inter. -Long 


0.55 


3x10-5 


Inter. 


4412 


2 546 


5 586 


Short-Long 


0.69 


< 10-6 


Long 


2 589 


1038 


7673 



Table 4: Left part: Results of the K-S test applied on the peak-count rates F for the RHESSI 
database. The K-S distance D and the K-S significance P are mentioned. Right part: The 
means, medians, and standard deviations of the peak-count rates are listed. 
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Table 5. Monte Carlo double-check of results from the statistical tests. 



Tests 


Inter.-Short 


Inter .-Long 


Short-Long 


Lags 


8 556 (85.6%) 


6 938 (69.4%) 


(0.0%) 


Norm, lags 


9 936 (99.4%) 


8 862 (88.6%) 


1458 (14.6%) 


Peak rates 


47 (0.5%) 


(0.0%) 


(0.0%) 



Note. — The number of cases out of 10000 MC cycles (and 
their percentages) arc noted for A-D (lags and norm, lags) 
and K-S (peak rates) probability values exceeding 5 % for tests 
done on spectral lags, normalized lags, peak-count rates, and 
on different pairs of GRB groups. 



Inter.-short 
inter.-long 
shorl-long 



4000 6000 
sorted order 



inter.-short 
Inter.-long 
short- long 



4000 6000 
sorted order 



Fig. 6. — Left panel: The A-D probabihties of the tests apphed on the samples of lags (left 
panel) and normalized lags (right panel) obtained from 10 000 MC cycles for different GRB 
groups. The horizontal solid hue denotes the 5 % threshold. 
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Fig. 7. — The K-S probabilities of the tests appUed on the samples of peak-count rates 
obtained from 10 000 MC cycles for different GRB groups. 
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Fig. 8. — Cumulative distributions of the spectral lags (left panel) and peak-fluxes (right 
panel) for the three BATSE GRB groups are shown. 



Groups 


A-D P 


Group 


Mean L 


Median L 


a 




(%) 




(ms) 


(ms) 


(ms) 


inter.-short 


51.3 


short 


177.1 


72.0 


454.6 


inter.-long 


3.8 


inter. 


207.5 


60.0 


464.2 


short-long 


9.7 


long 


390.7 


94.0 


848.2 



Table 6: Left part: Results from the A-D test of the equality of the spectral lag distributions 
for the BATSE GRBs. P denotes the P- value of the test. Right part: The mean, median 
and standard deviations a of the lags. 
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Groups 


D 


K-S P 


Group 


Mean 


Median 


a 






(%) 




F 


F 




inter.-short 


0.30 


< 10-6 


short 


4.00 


2.15 


6.80 


inter. -long 


0.13 


1.0 


inter. 


3.15 


1.29 


4.91 


short-long 


0.21 


< 10-6 


long 


4.09 


1.51 


10.31 



Table 7: Left part: Results of the K-S test applied on the peak-fluxes F (ph.cm-^.s-^) of 
the BATSE GRBs. The shortcuts have the same meaning as in Tab. HI Right part: The 
means, medians, and standard deviations of the peak-fluxes are listed. 




2468 2468 
number of components number of components 



Fig. 9. — Left panel: Bayesian information criterion (BIG) values for different models as a 
function of the number of bivariate Gaussian components. The higher the BIG value, the 
more probable the model. The most probable model is EVI with two components. The data 
sample consists of two variables: Tgo durations and hardness ratios. Right panel: BIG values 
for different models plotted against the number of components. The most probable model 
is EEE with three components. The data sample consists of three variables: Tgg durations, 
hardness ratios H, and peak-count rates F. 
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Fig. 10. — A scatter plot of 427 bursts, with measured T90, H, and F assigned into three 
groups by the EEE model. 



Table 8. A summary of the results from the model-based clustering. 





Model 


k 


BIG 


ABIC 


ABIC 


ABIC 


ABIC 


Evidence 










k=l 


k=2 


k=3 


k=4 




2 par. 


EVI 


2 


-681.5 


> 10 


X 


> 10 




very strong 


3 par. 


EEE 


3 


-1156.6 




> 10 


X 


> 10 


very strong 


4 par. 


VVV 


2 


-1768.4 


> 10 


X 


> 10 




very strong 



Note. — The results for model-based clustering applied on 2, 3, and 4 parameters is 
presented. The values of BIC for the best fitted models with k components are listed, 
as well as the differences to the models with other number of components. 
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Table 9. The spectral lags and peak-cmmt rates of the RHESSI GRBs. 



GRB=- Groupb L (ms)'= F (s-i)d ap {s-'^T 
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289.6 

4.36.2 
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Table 9 — Continued 



GRB* Groupb L (ms)<= F (s-i)d ap (s-i)« 
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Table 9 — Continued 



GRB^ Group'' L (ms)'^ F (s''^)'^ ap (s"!)" 
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Tablc 9 — Contimicd 



031118 


1 




4669.8 


592.7 


031120 


3 


1147.6 tweL 


991.2 


57.3 


031127 


3 




634.4 


94.5 


031130 


3 




1354.9 


155.0 


031214 


3 




609.8 


159.3 


031218 
031219 
031226 


1 
3 
3 


,3 - +137.0 

334.0 1589-0 


4732.5 
3262.1 
640.4 


1067.7 
439.7 
72.5 


031226 


3 




612.8 


86.8 


040102 


3 




1160.1 


139.1 


040108 


3 




395.7 


70.3 


040113 


3 




1183.6 


193.7 


040115 


3 




473.1 


84.0 


0401z5 


3 




528.6 


111.7 


040z0oA 


3 




497.4 


152.7 


U4UzUO-t5 


Q 
O 




oio.l) 


10.0 


ri/inoriT 
U4UzU ( 


Q 
O 


no c +91.0 
23.6 4 


00 ( o.y 


110. y 


n /I no T 1 


3 




751.8 


181.5 


040215 


3 




391.7 


46.9 


040220 


3 


1041 1 +1167.5 
LU11..L _1171.1 


1714.0 


178.9 


040225A 


3 




576.5 


104.0 


040225B 


3 




642.5 


90.0 


n /I nono A 

040o02A 


q 
o 

3 


1Q 9 +39.5 
±3.^ -36.0 


11/1 CQ C 

712.0 


zDl.D 

98.2 


040302B 


3 


ini Q +47.2 
101.8 I32.4 


13904.2 


283.4 


040303 


3 




403.8 


86.2 


040312 


1 


9 Y +27.4 
-14.8 


7197.9 


690.1 


040316 
040323 


3 
3 


1 fil Q +316.9 


4356.0 
650.8 


452.4 
174.9 


040324 


1 


9 7 +4.5 


16984.5 


1112.0 


040327 


3 




405.8 


56.2 


040329 


2 


fj +9.2 
•J-" -9.6 


19974.8 


849.1 


040330 


3 




584.5 


98.3 


040404 


3 




1096.6 


183.6 


040413 


1 




5163.3 


594.7 


040414 


3 




992.3 


59.5 


040421 
040423 


3 
3 


98.6 trti 


7197.3 
663.4 


220.2 
111.5 


040425 
040427 


3 

3 


155.3 I126.2 


7424.3 

674.8 


843.8 
112.1 


040429 


3 




582.8 


85.4 


040502A 


3 




3551.8 


310.5 


040502B 


3 




602.9 


48.9 


040506 


3 




807.0 


61.2 


040508 


3 




242.0 


64.4 


040510 


3 




1097.3 


123.3 



-31 - 



Tabic 9 — Contimicd 



GHB" (iioui)'' L (ms)' (s ay (s \)'' 



040513 


3 




218.9 


49.2 


040526 


3 




499.1 


96.2 


040528 


3 


_ +1067.9 
loJ-O _ii9i.3 


1827.0 


119.2 


040531 


3 




1024.5 


75.3 


040601 


3 




401.3 


93.9 


040603A 


3 




726.3 


138.6 


040603B 


3 




161.9 


26.5 


040oOoA 


3 




800.3 


224.4 


040d05B 


1 


n +46.0 
26-0 Ibo.8 


8657.1 


1243.4 


040605C 


3 


1271.6 


87.0 


U4UDli 
040619 


3 
3 


ont;'7 n +2983.7 
-2267.0 _ 1345.5 


740.1 
3308.1 


87.8 
480.5 


040701 


3 


1 nc p +246.2 
100. U _3i8.9 


1128.0 


103.9 


040719 


3 




1445.1 


177.0 


040723 


3 




4168.6 


697.7 


040731 
040803 


3 
3 


lino +244.6 
110-2 I2B7.4 


1803.4 
267.3 


91.5 
39.1 


040810 


3 




1654.2 


109.7 


040818 
040822 


3 
2 


qq Q +97.4 

^•^■^ -153.4 


2681.5 
2226.7 


225.6 
374.5 


040824 


3 




282.2 


45.8 


040921 


1 




2619.0 


490.6 


040925 
040926 
041003 


3 
3 
3 


1 nr- c\ +2908.9 

A 1 +92 7 
^•^ -83.7 


1263.8 
5855.8 
735.4 


167.1 
207.3 
182.2 


041006 


3 




1307.7 


254.1 


U41UU7 

041009 


2 

3 


ini n +119.5 
101.0 Ii23.3 


3301.3 

1382.6 


238.6 

111.5 




1 


n. 7 +6.3 
-8.0 


looo.o 


DID. Z 


041012 


3 




276.6 


44.5 


041013A 


3 




530.3 


39.7 


041013B 


1 




3412.2 


421.8 


041015 


3 




1630.8 


189.7 


041016 


3 




509.0 


78.5 


041018 
041101 


3 
3 


fi774 4 +6766.0 
D( 11.1 _7424.5 


717.6 
1036.3 


120.4 
167.5 


041102 


3 




1437.3 


188.3 


041107 


3 




883.0 


83.3 


041116 


3 




269.7 


52.1 


041117 


3 




1165.2 


106.5 


041120 


3 




1167.1 


123.6 


041125 


3 


-7.8 tt\-l 


9489.2 


168.7 


041202 


3 


76.3 t'^,f 


5824.1 


198.0 


041211A 


3 


Kf< ^ +667.3 
oo.o _825.0 


1066.9 


83.9 


041211B 


3 




5126.6 


629.1 


041211C 


3 


-6 +^^'-' 
-12.9 


11832.5 


315.0 



-32 - 



Tabic 9 — Contimicd 



(;HI!-' Gnmp'' L (ms)' I' ^ at (s M 



041213 


1 




6872.1 


669.2 


041218 


3 




364.9 


54.7 


041219 


3 




560.3 


76.9 


041223 
041224 


3 
3 


1 f. +208.1 
-l.b _2ii.7 


1567.8 
285.4 


91.6 
73.2 




2 


Ol fi +187.2 
^l.D _i77.7 


2211.1 


206.3 


U5Ulz4 


3 




815.2 


162.3 


05012D 


3 




1477.0 


79.6 


050203 


3 




2870.6 


234.3 


050213 

050214 


3 

3 


AA 9 +197.8 
^^■^ -194.9 


2174.0 

390.4 


101.5 

60.3 


U5U21D 
050219 
050311 


1 
3 
3 


f- n +99.4 
6-0 -76.6 
/lir» o +321.9 


4014. 1 

4760.8 
465.1 


bob. 4 
259.7 
88.4 


050312 


1 


4 4 +44.5 
-31.0 


4439.8 


419.4 


050314 


3 




2101.0 


151.0 


050320 


3 




838.7 


114.5 


050321 


3 




880.5 


128.4 


05032b 


3 




3501.9 


403.9 


050328 

050404 
050409 
050411 


1 

3 
2 

3 


,„ c +79.4 

-43.5 1 

r»n T +35.6 

^^•^ -39.0 
-1 Y +16.9 
-14.1 


12706.3 

6336.6 
18969.8 
813.9 


1069.6 

215.9 
1936.7 
125.9 


050412 


3 




1455.7 


140.8 


050429 


3 




2018.6 


145.8 


050430 


3 




986.8 


139.4 


050501 


3 




1697.9 


308.8 


050502 


2 




1132.2 


235.6 


050509 


3 


^ -I-3Q7 4 


2282.4 


145.7 


050516 


3 




259.7 


71.5 


050525A 


3 




4853.5 


259.6 


050525B 

050528 


3 

3 


o c +60.0 
-50.3 


7311.6 

421.3 


214.3 

48.7 


050530 


3 




1963.1 


273.3 


050531 
050614 


3 
3 


56.6 lig3 g 


6082.8 
377.5 


197.2 
58.6 


050701 


3 




1173.5 


166.1 


050702 


2 




951.9 


187.8 


050703 


3 




2734.2 


305.4 


050706 


3 




1620.7 


168.6 


050713A 


3 


QQQ Q +1804.9 
-ZOO.O _2189.4 


1450.6 


244.4 


050713B 


3 




431.1 


69.1 


050715 


3 




2053.4 


198.9 


050717 


3 


110 3 +151.6 
J.1U.O _220.0 


2456.3 


143.1 


050726 


3 


810 5 +"51.7 
oiu.o _1334.3 


1480.4 


90.7 


050729 


3 




873.7 


177.5 



-33- 



Table 9 — Continued 



GRB* Group'' L (ms)<= F (s-i)d ap 



UoUoUz 


3 




498.5 


97.6 


050805 


2 


(■ Q +30.0 
-20.0 


3355.3 


334.6 


uououy 


Q 
O 




4oyc5.o 


oyy.o 


OoOolo 


3 




810.7 


110.4 


OoOol4 


1 




0499.0 


1658.3 


UoUolT 


3 




895.5 


118.2 


OoOozU 


3 




681.3 


101.9 




1 


1 o +3.2 
1-2 -4.0 


7474.8 


930.5 


050825 
Oooyoi 


2 
3 


on a +132.9 
-82.6 3 


4307.4 
631.9 


585.3 
149.5 


U5U92d 


3 




1872.2 


395.3 


051009 


3 




475.8 


66.9 


051012 
051021 
051031 


3 
3 
3 


/IRO c; +443.7 
462.0 _207.9 
Q7 S +716.9 
-y-o -496.2 


/80.0 
2117.7 

678.7 


331.8 
221.7 
45.4 


051101 


3 




921.2 


258.0 


051103 


1 


n R +2.4 

0.6 I2.8 


135199. D 


10873.8 


051109 


3 




882.0 


92.0 


051111 


3 




516.1 


102.1 


051117 


3 


OO/I « +2133.4 
-^z*.o _33i8.9 


942.4 


133.9 


051119 


3 




7 43 A 


133.7 




Q 
O 








UollZ4r> 


3 




2715.0 


152.8 


n cc 1 on 1 A 


3 




211.7 


37.0 


051201B 


3 




1000.9 


170.4 


UoizU / 


Q 
O 






1 01 


051211 


3 


CQO q +836.1 
dsz.o _gi4.9 


1573.5 


69.6 


051217 


3 




458.7 


81.7 


051220A 
051220B 


3 
3 


39.0 


19918.0 
321.2 


658.7 
73.1 


051221 


1 


0.0 


15280.2 


1298.8 


051222 


3 




334.4 


78.5 


060101 
060110 


3 
3 


599.3 


1647.1 
445.6 


94.2 
107.8 


060111 
060117 


3 
3 


-2346.6 l?03i.i 


1079.4 
1300.6 


46.7 
113.7 


060121A 


3 




2980.4 


277.6 


060121B 


3 




343.0 


46.7 


060123 


3 


IS 4 +624.7 
-529.7 


3965.2 


184.2 


060124 


3 




713.2 


76.3 


060130 


3 




2285.1 


376.4 


060203 
060217 


1 
3 


90 t; +46.5 


6089.8 
2432.4 


831.6 
416.0 


060224 

060228 


3 

3 


q4A n +3150.4 

Tno n +1896.9 
-(UZ.U _i770.9 


851.0 
717.4 


109.2 

89.8 


060303 


1 


91 9 +46.3 
^1-^ -53.5 


9343.4 


999.5 



-34- 



Table 9 — Continued 
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A. Uncertainties of Decimated Data 



This part describes the one sigma uncertainties of the bin counts for the decimated 



RHESSI data (ICurtis et al. II2OO2I : ISmith et al. Il2002f ). 



A.l. Full decimation 

Calculation of one sigma uncertainty ac^^ for the bin counts Cdc of fully decimated data 
and then corrected for this decimation is the following. For corrected bin counts Cdc it holds: 

Cdc = /d.Cd, (Al) 

where /d is the decimation factor (weight), usually equal to 4 or 6 for the RHESSI data, and 
Cd is the number of counts in a bin of the decimated signal. If we assume that counts in a 
bin follow Poisson statistics then: 

o"Cd = /d^Cd = a/ /d-Cdc, (A2) 
where ac^ is the dispersion of the Cd. 



da 



dc 



A. 2. Partial decimation 

Now consider the situation in which counts in a bin are only partially decimated, i.e. 
they consists of the non-decimated Ci and the decimated signal C2,d- This situation may 
happen when we sum counts over the energy band [Ei, E2], Ei < Eq < E2, and only the 
counts below the energy Eq are decimated. Then the corrected signal Cdc is equal to: 



Cdc — Ci + C2_dc — Ci + /d-C2,d; 



(A3) 
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Tablc 9 — Continued 
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674.6 


120.5 



'^RHESSI GRB number. 

''The assignment to the GRB group; 1 - short, 
2 - intermediate, 3 - long. 

'^Spectral lags were calculated from the difference of the 
count light-curves at the energy intervals 400 — 1500 keV and 
25 — 120 keV. The errors compose of the 95 % CL statistical 
uncertainty and the light-curve's time resolution. 

Peak-count rates derived in the band 25 — 1500 keV. 

"One sigma statistical uncertainties of the peak-count rates. 
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where /d is again the decimation factor, and C2,dc is corrected part of the signal that was 
decimated. One sigma uncertainty ac^^ is then equal to: 

/ / 3C \ ^ f dC \ ^ 

•"^^ = V [m ^ (s^J = 

= sj-'h + /d-<,. = + A'-C2,d = VC, + /d.C2.d„, (A4) 

where o-f^ = ^/Ci is the dispersion of the non-decimated part of the bin counts and 0"C2d ~ 
\/C2,A is the dispersion of the decimated part of the bin counts. 
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